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Abstract. A plasma velocity shear layer and/or a tangential flow discontinuity provide conditions 
allowing for energetic particle acceleration. We review such acceleration processes acting both in non- 
relativistic and in relativistic flows. In heliospheric conditions shear layers can provide particles with 
energies compatible with the observed values (from several keV up to MeV), while in relativistic 
extragalactic jets proton energies even in excess of 10 19 eV can be obtained. Application of the 
discussed theory to particular astrophysical objects is severely limited by inadequate knowledge of 
local physical conditions. 



1 Introduction 

The first order Fermi acceleration in shock waves and the second order acceleration in 
turbulent MHD media are widely considered as main sources of cosmic ray particles in 
astrophysical conditions. In the present paper we consider an alternative mechanism 
involving particle acceleration at velocity shear layers formed in non-uniform plasma 
flows, e.g. in a magnetosheath enveloping the Earth magnetosphere or at the interface 
between the relativistic jet and its ambient medium. Till now the considered compli- 
cated physical phenomenon was only occasionally discussed in the literature. The 
process was introduced into consideration by Berezhko and collaborators in a series 
of papers in early eighties (cf. Berezhko 1981, 1982a, b; Berezhko & Krymsky 1983; 
Bczrodnykh et al. 1984a, b, 1987; summarised in a review by Berezhko 1990). Much 
later an independent discussion of such processes acting in non-relativistic shear lay- 
ers was presented by Earl et al. (1988), Jokipii et al. (1989) and Jokipii & Morrill 
(1990), and for relativistic tangential flow discontinuities by Ostrowski (1990; cf. also 
Berezhko 1990). A discussion of possible cosmic ray acceleration in mildly relativistic 
jets up to ultra-high energies and consequences of acting such acceleration processes 
in ultra-relativistic ('mili-arc-second') jets were considered in recent papers by Os- 
trowski (1998a, b; 1999). Below, we will shortly discuss the main results obtained by 
the above authors. 
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2 Particle acceleration in a shear layer 

A high energy particle scattered after crossing a shear flow layer can gain or loose 
energy. It is due to a respective velocity difference of the final scattering centre rest 
frame with respect to the particle starting point, 

AU^^Ax, (1) 
ax 

where we consider a 1-D situation with the flow velocity U directed along the z- 
axis, and the velocity gradient along the x-axis of the reference frame. In absence 
of magnetic field Ax = v x At is a free path along the x-axis (v — [v x ,v y ,v z ] is the 
particle velocity). Let us assume for a while the scattering centres to be static with 
respect to the local plasma rest frame. Then, in the scattering centre rest frame the 
particle momentum changes with respect to the one in the starting point plasma rest 
frame at 

Ap=^. (2) 

V 

For a mean AU <C v the full process can be described as the momentum diffusion 
with the diffusion coefficient 

i/(A P r\ v 2 (du\\ 

where the second equality comes from averaging over an isotropic particle distribution, 
t = (At) is the mean scattering time and the term (dU/dx) 2 is the shear scalar in 
the considered simple flow pattern. 

In the presence of magnetic field the mean particle shift in the x direction can 
be much smaller than v At. Then the introduced r parameter equals the ratio of the 
particle mean free path (shift) along the x-axis, X x , to the respective mean parti- 
cle velocity (v x ), t — X x / (v x ). More exactly the particle energy change and the r 
parameter in Eq. 3 has to be derived by averaging over actual particle trajectories. 

If the parameter t scales with the particle momentum as r cx p n , then the ac- 
celeration process acting within the shear layer produces the high energy asymptotic 
phase-space distribution (Berezhko 1982) 

/(p)«p-( 3+ ">. (4) 

One should note that the considered acceleration process in a shear layer plays a sub- 
stantial role if the mean plasma velocity difference at successive scatterings (Eq. 1) 
is larger than the turbulent velocities leading to the ordinary second-order Fermi 
acceleration (cf. Eq. 10 below). 

3 Cosmic ray viscous acceleration in the Heliosphere 



Because of insufficient information about the turbulent shear flow patterns within the 
astrophysical shear layers it is not possible to give a firm evaluation of the viscous 
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acceleration rates. However numerous measurements show enhancements of energetic 
particle populations in the Hcliosphere, where the Solar Wind forms shear flows. In 
these cases one can estimate the highest energies for accelerated particles by the 
viscous mechanism and compare its to the measured ones. Within the Heliosphere 
such estimates were provided for the observed shear flow sites (cf. Berezhko 1990, 
Jokipii & Morfill 1990), including the Earth and the Jupiter sheared magnetosheath, 
interplanetary magnetic field sector boundaries, boundaries of the high speed Solar 
Wind streams. In general the evaluated energies are within the observed ranges. 

4 Particle acceleration at relativistic shear layers 

The relativistic shear layers occur in a number of objects in space, including galac- 
tic and extragalactic relativistic jets and accretion discs near black holes. Below we 
consider the jet side boundary layer as an example of the relativistic shear flow. 

For particles with sufficiently high energies the transition layer between the jet 
and the ambient medium can be approximated as a surface of a discontinuous ve- 
locity change, a tangential discontinuity ('td'). It becomes an efficient cosmic ray 
acceleration site provided the considered velocity difference U is relativistic and the 
sufficient amount of turbulence is present in the medium. The situation with highly 
relativistic jet (T = (1 — [/ 2 )~ 1//2 ^> 1) was not quantitatively discussed till now 
and, thus, our present discussion is mostly based on the results derived for mildly 
relativistic flows by Ostrowski (1990, 1998a). 

4.1 Energy gains 

Any high energy particle crossing the jet boundary changes its energy, E, according 
to the respective Lorentz transformation. It can gain or loose energy. In the case of 
uniform magnetic field the successive transformation at the next boundary crossing 
changes the particle energy back to its original value. However, in the presence of 
perturbations there is a positive mean energy change: 



The numerical factor tje increases with the growing magnetic field perturbations and 
slowly decreases for increasing T. For mildly relativistic flows, in the strong scattering 
limit particle simulations give values of t]e as substantial fractions of unity (Ostrowski 
1990). For large T we assume the following scaling 



where rj = f](T = 2). In general 770 depends also on particle energy. During the accel- 
eration process, particle scattering is accompanied with the jet's momentum transfer 
into the medium surrounding it. On average, a single particle with the momentum p 
transports across the jet's boundary the following amount of momentum: 



(AE)= m (T-l)E. 



(5) 
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(6) 



(Ap) = (A Pz ) = Vp (T - 1) Up, 



(7) 
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where the z-axis of the reference frame is chosen along the flow velocity. The numerical 
factor ri p s=s rj E and there acts a drag force per unit surface of the jet boundary and the 
opposite force at the medium along the jet, of the magnitude order of the accelerated 
particles' energy density. Independent of the exact value of r/E, the acceleration process 
can proceed very fast due to the fact that average particle is not able to diffuse - 
between the successive energizations - far from the accelerating interface. One should 
remember that in the case of shear layer or tangential discontinuity acceleration and, 
contrary to the shock waves, there is no particle advection off the 'accelerating layer'. 
Of course, particles are carried along the jet with the mean velocity of order U/2 and, 
for efficient acceleration, the distance travelled this way must be shorter than the jet 
breaking length. 

The simulations (Ostrowski 1990) show that the discussed acceleration process 
can be quite rapid, with the time scale given in the ambient medium rest frame as 



where r g is a characteristic value of the particle gyroradius. For efficient scattering the 
numerical factor a can be as small as ~ 10 (Ostrowski 1990). One may note that the 
applied diffusion model involves particles with infinite diffusive trajectories between 
the successive interactions with the discontinuity. However, quite flat spectra, nearly 
coincident with the stationary spectrum (cf. Fig. 1), are generated in short time scales 
given by Eq. 8 and these distributions are considered in the present discussion. For 
the mean magnetic field B g given in the Gauss units and the particle (proton) energy 
EecV given in EeV (1 EeV = 10 18 eV) the time scale (8) reads as 



For low energy cosmic ray particles the velocity transition zone at the boundary is 
a finite-width turbulent shear layer. We do not know of any attempt in the literature 
to describe the internal structure of such layer on the microscopic level (cf. Aloy 
et al. 1999, Hcnrikscn, at this meeting). Therefore, we limit the discussion of the 
acceleration process within such a layer to quantitative considerations only. From 
rather weak radiation and the observed effective collimation of jets in the powerful 
FR II radio sources one can conclude, that interaction of a presumably relativistic jet 
with the ambient medium is relatively weak. Thus the turbulent boundary layer must 
be relatively thin, with thickness denoted with D. Within it two acceleration processes 
energise low energy - the ones with the mean radial free path A <C D - particles. The 
first one, discussed in section 2 above, is connected with the velocity shear and is 
called 'cosmic ray viscosity'. The second one is the ordinary Fermi process in the 
turbulent medium. The acceleration time scales can not be evaluated with accuracy 
for these processes, but - for particles residing within the considered layer - we can 
give an acceleration time scale estimate 




(8) 



Ttd ~ IO-^eovSg 1 [s]. 



(9) 



2-nr, 



(10) 



m = 



c ^ + (C/A) 2 ' 
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Fig. 1. Spectra N(t, p, x = 0) = dn/d log p of accelerated particles at the jet boundary in a sequence 
of times ti < t2 < ■ ■ ■ < tj (each tj = 10ti— i), for continues particle injection at small momentum 
P0 < Pc- In the upper plot the spectra generated in the turbulent shear layer are given, while, in 
the two lower ones, the tangential discontinuity generated spectra arc presented: the results for the 
particle mean free path X oc p are given in the upper panel, and for A = const in the lower one. 
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where V is the turbulence velocity. One expects that the first term in the denominator 
can dominate at low particle energies, while the second one for larger energies, with 
Tu approaching the value given in Eq. 8 for A ~ D. If the second-order Fermi accelera- 
tion dominates, A < D(V/U), the time scale (10) reads as r n ~ 10 8 E TeV B^V^ 2 [s], 
where V3 is the turbulence velocity in units of 3000 km/s. Depending on the choice 
of parameters this scale can be comparable or longer than the expansion and internal 
evolution scales for relativistic jets. In order to efficiently create high energy particles 
for the further acceleration by the viscous process and the tangential discontinuity 
acceleration one have to assume that the turbulent layer includes high velocity tur- 
bulence, with V3 reaching values substantially larger than 1, or other high energy 
particles sources are present. For the following discussion we will assume that such 
effective pre-acceleration takes place, but the validity of this assumption can be es- 
timated only a posteriori from comparison of our conclusions with the observational 
data. 

4.2 Energy losses 

To estimate the upper energy limit for accelerated particles, at first one should com- 
pare the time scale for energy losses due to radiation and inelastic collisions to the ac- 
celeration time scale. The discussion of possible loss processes is presented by Rachen 
& Biermann (1993). The derived loss time scale for utra-high energy protons can be 
written in the form 

T loss ~5-10 9 B-* (1 + Xa)- 1 E£ w [s], (11) 

where a is the ratio of the energy density of the ambient photon field relative to 
that of the magnetic field and X is a quantity for the relative strength of p"f in- 
teractions compared to synchrotron radiation. For cosmic ray protons the acceler- 
ation dominates over the losses (Eqs. 9, 11) up to the maximum energy -EecV ~ 
2-10 2 a- 1 [£ G (l + A:a)r 1/2 - 

4.3 Spectra of accelerated particles 

The acceleration process acting at the tangential discontinuity of the velocity field 
leads to the flat energy spectrum and the spatial distribution expected to increase 
their extension with particle energy. Below, for illustration, we propose two simple 
acceleration and diffusion models describing these features. 

A turbulent shear layer 

At first we consider 'low energy' particles wandering in an extended turbulent shear 
layer, with the particle mean free path A oc p. With the assumed conditions the 
mean time required for increasing particle energy on a small constant fraction is 
proportional to the energy itself, and the mean rate of particle energy gain is constant, 
(^)gain = const. Let us take a simple expression for the synchrotron energy loss, 
(E)i oss oc p 2 , to represent any real process acting near the discontinuity. One may 
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note that the jet radius and the escape boundary distance provide energy scales to 
the process. Another scale for particle momentum, p c , is provided as the one for equal 
losses and gains, (E) galn = (E) loss . As a result, a divergence from the power-law and 
a cut-off have to occur at high energies in the spectrum. 

We use a simple Monte Carlo simulations to model the acceleration process for 
a continuous particle injection, uniform within the considered layer. The diffusion co- 
efficient k±_ is taken to be proportional to particle momentum, but independent of the 
spatial position x. We neglected particle escape through the shear layer side bound- 
aries and we assumed §^ = . For the escape term we simply assume a characteristic 
escape momentum p max . In Fig. 1 we use p c as the unit for particle momentum, 
so it defines also a cut-off for p c < p max . At small momenta the spectrum has a 
power-law form - in our model the averaged over angles f(t,p) oc p (cf. Eq. 4) 
- with a cut-off momentum growing with time. However, in long time scales, when 
particles reach momenta close to p c , losses lead to spectrum flattening and pilling up 
particles at p close to p c . Then, a low energy part of the spectrum does not change 
any more and only a narrow spike at p p c grows with time. Let us also note that 
in the case of efficient particle escape, i.e. when p max < p c , the resulting spectrum 
would be similar to the presented by Ostrowski (1998a) short time spectrum with 
a cut-off at « p max . 

Tangential discontinuity acceleration 

An illustration of the acceleration process at the tangential discontinuity have to take 
into account a spatially discrete nature of the acceleration process. Here, particles are 
assumed to wander subject to radiative losses outside the discontinuity, with the mean 
free path oc p and the loss rate oc p 2 . At each crossing of the discontinuity a particle 
is assumed to gain a constant fraction A of momentum (cf. Eqs. 5, 6), p' = (1 + A)p , 
and, due to losses, during each free time At its momentum decreases from p m to p 
according to 

-- — = const • At. (12) 

V Pin 

The time dependent energy spectra obtained within this model are presented in lower 
panel in Fig. 1, where we choose units in a way to put the constant in Eq. 12 equal 
to one and the particle mean free paths are equal in two considered models at p = p c . 
Comparison of the results in two models allows to evaluate the modification of the 
acceleration process by changing the momentum dependence of the particle diffusion 
coefficient. For slowly varying diffusion coefficient (the 'A = const' model) high energy 
particles which diffuse far away off the discontinuity and loose there much of their 
energy still have a chance to diffuse back to be accelerated at the discontinuity. In 
the model with k quickly growing with particle energy (the 'A = C ■ p' model) such 
distant particles will decrease their mobility in a degree sufficient to break, or at least 
to limit their further acceleration. One should note that in both models the spectrum 
inclination at low energies is the same (here the particle density n(p) oc p~ 2 ). 
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5 Final remarks 

Shear layers occurring in astrophysical plasma flows are able to accelerate cosmic 
ray particles in the so called viscous acceleration process. Depending on conditions 
the process can be described as the particle momentum diffusion or the tangential 
discontinuity acceleration. The last one can be very efficient in relativistic flows. In 
particular, in jets in active galactic nuclei the cosmic ray protons can reach energies 
in excess of 10 18 eV. 

The generated cosmic ray populations can influence the shear layer flow through 
viscous and/or dynamical forces (cf. Arav & Begelman 1992, Ostrowski 1999). In 
relativistic jets the so called cosmic ray cocoon can be formed leading to a number 
of observational effects. The essential problem with application and verification of 
the presented theory is insufficient information about the local parameters of the 
considered shear layers. One may note several recent observational papers showing 
effects which could be ascribed to, or are at least compatible with the acceleration 
process acting at jet boundary layer (Attridge et al. 1999, Scarpa et al. 1999, Perlman 
et al. 1999). 
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